Collisions between solitary waves of three-dimensional Bose-Einstein condensates 
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We study bright solitary waves of three dimensional trapped Bose-Einstein condensates and their 
collisions. For a single solitary wave, in addition to an upper critical number, we also find a lower 
cut-off, below which no stable state can be found. Collisions between solitary waves can be elastic, 
inelastic with either reduced or increased outgoing speed, or completely unstable due to a collapse 
instability. A 7r-phase difference between the waves promotes elastic collisions, and gives excellent 
agreement with recent experimental results over long timescales. 
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Bright solitary waves (BSWs) of atomic Bose-Einstein 
condensates (BECs) have been observed in the form of 
a single BSW [| and groups of interacting BSWs 0, Q 
for attractive atomic interactions, and as gap BSWs in 
eriodically-confmed BECs with repulsive interactions 
In the former case, the attractive interactions which 
support the BSW against dispersion also induce a col- 
lapse instability when the number of atoms exceeds a crit- 
ical value Nq 0-E3]- making attractive BECs problematic 
to form. Consequently, 'soliton' experiments first form a 
stable BEC with repulsive atomic interactions and switch 
to attractive interactions by means of the molecular Fesh- 
bach resonance 0. For N < Nc, a single BSW forms Q, 
while for N > Nc a 'bosenova'-type collapse instability 
|E B 13 is induced with excess atoms being ejected from 
the condensate. The 'remnant' condensate is observed 
to contain BSWs 0, Q which oscillate in the trap and 
interact with each other. It has been inferred Q] that 
the BSWs repel each other due to 7r-phase differences 
between them, thought to be formed by modulational in- 
stability of the condensate 0, IToL Hlj . With matter- wave 
'solitons' being potential candidates for applications such 
as interferometry and quantum information processing, a 
thorough understanding of their properties and in partic- 
ular their interactions is of key importance. 

Theoretically, bright solitons are stable solutions of the 
ID nonlinear Schrodinger equation with a sech-shaped 
wave profile [12. In isolation their form is preserved by 
a balance between dispersion and attractive nonlinearity, 
and they can have any population. Another key property 
is that they emerge from collisions with unchanged form. 
However, even in ID, the exact collisional form depends 
on the relative phase between the incoming solitons A<p 
(l3j| : for A<p — the waves can overlap completely (at- 
tractive solitonic interaction), while for Acf> = ir no over- 
lap can occur and the solitons appear to 'bounce' (re- 
pulsive solitonic interaction). BSWs, the 3D analog of 
the bright soliton, may have considerably different prop- 
erties: the existence of multiple dimensions leads to the 
collapse instability 5,9] and thermal dissipation [l4| . In- 
deed, two colliding BSWs can form a high-density state 
unstable to collapse, with relative phase difference play- 
ing a key role [Tl], Ha] • Whereas the experiments of 0, Q 



use highly-elongated 'quasi-lD' geometries, Cornish et al. 
employ a very weakly-elongated geometry, and, some- 
what surprisingly, observe well-behaved BSW oscillations 
for very long times with negligible damping. This raises 
questions over how 'solitonic' such 3D states are. 

In 3D an untrapped attractive BEC is always unstable 
to collapse 01- Addition of confinement can stabilise it 
up to a critical density and number of atoms Nq [H, UtI 
ITsl IT^ ] . By definition a BSW should be capable of axial 
self-trapping. In the absence of external axial trapping, 
radial confinement is necessary to make Nq finite, and 
also to enable approximation to a ID soliton. 

In this work we simulate BSWs of 3D BECs and their 
collisions. We observe rich and complex dynamics, which 
extend well beyond the ID soliton. Below a lower critical 
atom number we show that stable BSW solutions do not 
exist. Collisions can be elastic, inelastic with increased or 
decreased speed, or completely destructive. We highlight 
the effect of relative phase, number of atoms and speed on 
the collisional stability. Furthermore we show that recent 
observations of Cornish et al. |3[ are consistent with two 
oscillating BSWs with 7r-phase difference, with the simu- 
lations giving excellent agreement with the experimental 
oscillation data, over long timescales (for over 20 oscilla- 
tions/3 seconds), despite neglecting higher-order effects 
such as thermal dissipation |14j. 

We employ the Gross-Pitaevskii equation (GPE) to de- 
scribe the condensate mean-field 'wavef unction' ^(r, i), 
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Here m is the atomic mass and g — Airh 2 a§ / 'm is the 
interaction coefficient, where as is the s-wave scattering 
length (as < in this study, corresponding to attractive 
atomic interactions). The harmonic confining potential 
is cylindrically-symmetric with z and r representing the 
axial and radial coordinates. We base our analysis on the 
85 Rb experiment with scattering length as = —0.6 nm 
and radial trap frequency w r — 2n x 17.5 Hz. 

We begin by considering isolated BSW solutions. To 
generate such states we propagate the GPE in imaginary 
time subject to the desired atom number Ns and exter- 
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FIG. 1: (a) Upper critical number Ac as a function of uj z 
(uj r = 2n x 17.5 Hz). Indicated are the lower critical number 
Nc for uj z — (dashed line) and the condensate regimes, 
(b) Radially-integrated axial density versus time for a BSW 
(N s = 1500) with lo z = 0. (c)-(d) Same as (b) but with 
initial density rescaled to (c) Ns = 2700 and (d) iVs = 1000. 
Dark/light regions represent high/low density. 



nal potential [2p|. Following previous studies 0, 0] 
we isolate the critical number for collapse Nq when the 
GPE no longer converges to a solution. Figure Qfa) 
shows the critical number Nc as a function of axial trap 
frequency uj z . For an axially- homogeneous system we 
find a range of values of N$ where solutions exist which 
are axially self-trapping (Fig.^b)), i.e. BSW solutions. 
For Ns > N c — (2450 ± 30) the wavepacket is unsta- 
ble to interaction-induced collapse, followed by explo- 
sive dynamics (Fig.^c)). In addition to the upper crit- 
ical number we also find a lower critical regime when 
N s < N c = (1280 ± 20). For N s < N c the wavepacket 
is no longer self-trapped axially but spreads over time 
(Fig.njd)), since the attractive interactions are no longer 
sufficient to prevent dispersion of the wave. Both the 
low- TVs and high-TVs instabilities do not occur for the ID 
soliton, which is robust to any population. 

Under axial confinement Nc decreases weakl y w ith ui z . 
The axial trapping increases the peak density |2l| . mak- 
ing the BSW more prone to the collapse instability. For 
lo z = 2tt x 6.8 Hz Q we find N c = (2305 ± 30), in 
excellent agreement with experimental observations of 
Nc ~ (2284 ±300) 6]. Alternative approaches have pre- 
viously over or under estimated Nc 0, 0] ■ 

Axial confinement provides a sufficient restoring force 
to axially trap the condensates and achieve stable solu- 
tions for A~ < N' c . However, such a state can no longer 
be classed as a BSW because it is not capable of axial 
self-trapping; it is then simply the ground state trapped 
attractive BEC. The regimes of the attractive system, i.e. 
instability to collapse, stable BSW, and trapped BEC, 
are indicated in Fig.^a). The states observed in Q ap- 
pear to have a population in the range N' c < N < Nc, 
consistent with being BSWs. 

We now consider collisions between two BSWs, each 
with Ns = 1500, in an axially-homogeneous system. In 
ID this collision is stable and elastic. Initially the BSWs 
are well-separated and given an axial momentum kick, 
ip = \ip\ expiimvtnZ / K) , such that they propagate towards 
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FIG. 2: (a)-(b) w Du t/win for a BSW collision with Acj) = 
(crosses) and Acj) = n (circles) in a iii 2 =0 {u r = 2tt x 
17.5 Hz) system as a function of (a) incoming speed Din (Ns = 
1500), and (b) atom number TVs (vin = 0.12 mms -1 ). The 
dashed line represents ^out — ^in- 

In the shaded regions BSW 
solutions are not accessible (light region: iVs < N^; dark 
region: iVs > Nc)- (c) Axial density (integrated over r) as 
a function of Vi n t for (i) A(f> = and Ui n = 0.03mms _1 , (ii) 
Acj> — and v- ln = 0.09 mms -1 , (iii) Acj> — and Vi n = 
0.24 mms -1 , and (iv) A(f> = 7r and v- lxl = 0.09 mms -1 . Each 
box represents the region [0,80 /J,m/v- m ] x [—40,40] /im. 



each other with speed Vi n . In isolation the BSWs prop- 
agate with constant speed and shape (see early times in 
Fig. GJc)). We also control the phase difference Acf> be- 
tween the BSWs by initially multiplying one BSW by 
cxp(zA(/>). We primarily consider the extreme cases of 
Acj) = and n. We consider speeds of the order of 
those observed experimentally assuming an ampli- 
tude zo ~ 10 /jm in the trap, harmonic motion suggests 
a maximum speed v — lu z z ~ 0.4 mms -1 . 

We parameterise the stability of a BSW collision in 
terms of the outgoing speed v out - Figure EJa) plots the 
ratio v ou t / Win as a function of incoming speed vi n for the 
cases of Acj) = (crosses) and Acj) = ir (circles). For 
Acj) = and low v in (v in <0.09 mm s -1 ), the BSWs do 
not emerge from the collision (Fig. EJc)(i)) 22]. The 
overlap of the BSWs forms an intermediate state with 
sufficiently high density that a destructive bosenova-like 
collapse instability is induced 0,0,0]. For Acj) = and 
intermediate i>i n (0.1 ^v- m <0.14 mms -1 ), BSWs emerge 
from the collision but with w out < V{ n . The density again 
reaches a critical level, but only a partial collapse in- 
stability is induced (Fig. E^c)(ii)). This leads to an ir- 
reversible transfer of energy from centre-of-mass kinetic 
energy to excite 3D collective modes of the BSWs. As the 
Vi n is increased further the collision becomes more elastic 
(v ut —> v in) and the shape excitations become negligible. 
These elastic dynamics, shown in Fig.|2Jc)(iii), illustrate 
how the BSWs pass through each other and form 3 fringes 
in a stable Acj) — collision. This limit is analogous to 
the corresponding ID soliton collision. 

In contrast, Acj) — it collisions are elastic throughout 
this range of w; n . Overlap of the BSWs at the centre of 
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FIG. 3: (a)-(b) Ratio « ou t/«in of a BSW collision with Acb = 
(crosses) and Acb = n (circles) in a ui z = system as a 
function of (a) «i n (for u r = 2n x 35 Hz and iVs = 1500) and 
(b) u r (for v in = 0.06 mms -1 and iV s = 1500). The dashed 
lines plot v ou t ~ Vin- Inset in (a): axial density (integrated 
over r) as a function of Vi n t for Acb = 0, Vi n = 0.06 mms -1 and 
uj r = 2-7TX35 Hz. The box represents the region [—40, 40] /imx 
[0,50 /im/nin]. In the shaded regions of (b) BSW solutions 
are not accessible (light region: Ns < Nq] dark region: iVs > 
N C ). 

mass is prevented and the BSWs 'bounce' (Fig.|5{c)(iv)). 
This restricts the increase of the peak density during the 
collision [2^, and stabilises against instability [Tl| . 

From the solutions generated in Fig. ^a) for lu z = 
we estimate the critical density for collapse to be nc ~ 
5 x 10 19 m -3 . All of the Ac/) = collisions in Fig. Ufa) 
exceed this density, while the Ac/> = ir collisions do not. 
The stabilisation of collisions with speed suggests that 
the timescale plays a key role. Indeed, in the elastic (in- 
clastic) regime v ln > 0.15 mms -1 (v- m is 0.15 mms -1 ) the 
time over which nc is exceeded becomes consistently less 
(greater) than around 5 ms, which is the characteristic 
collapse time observed experimentally 0. 

The stability is also sensitive to atom number Ns- 
Figure EJb) shows v ou t / v ln as a function of As (v- m = 
0.12 mms -1 ). Addition of atoms makes the BSW more 
prone to the collapse instability, destabilising the colli- 
sions with increasing TVs ■ For Ac/> = the collisions break 
down for low Ns {Ns ~ 1400) while the Acj> = it collision 
remains elastic up until just below Ac- 

For < Acb < 7r the degree of overlap during the col- 
lision varies between full overlap (Ac/) = 0) and no over- 
lap (Acb — n). Population transfer occurs between the 
BSWs, with the outgoing speeds becoming asymmetric 
to conserve momentum, complicating the dynamics. Al- 
though we concentrate here on Acf> — and tt we have 
tested their robustness by considering 10% deviations of 
each (Acb = O.lir and 0.97r), and observe that the regimes 
of stability are robust to such variations [24j . 

We now double the radial trap frequency to u> r = 
2ir x 35 Hz. The collisional stability as a function of 
incoming speed is shown in Fig. [2Ja) . For low U; n we 
again observe inelastic collisions, but this time v out > v- m 
(^out/win ^4.5). As the incoming speed is increased fur- 
ther the collision becomes more elastic (iwt — ► v- ul ), as 
seen previously in Fig. Ola). A low speed inelastic col- 
lision with i>out/vin ~ 3.3 is shown in Fig. |3Ja)(inset). 



During the collision the high density triggers a partial col- 
lapse instability. A small fraction of atoms (around 10% 
in this case) are ejected from the BSWs, releasing a sig- 
nificant amount of interaction energy which is converted 
into kinetic energy. For lower cu r this would tend to gen- 
erate a radial burst of highly-energetic atoms (bosenova) 
0. However for sufficiently high to r we propose that 
this radial burst is suppressed, with a significant fraction 
of this kinetic energy focussed axially into the BSWs, 
thereby increasing their speed. Again, for Acb = n the 
collisions remain elastic throughout this range of speeds. 

In Fig. EHb) we vary the radial trap frequency u> r 
while keeping Ns fixed. We observe BSW destruction 
for Acb = and low u r (3D bosenova instability), while 
for larger Lo r we typically see the ratio i> ut/^in grow. 
There is occasionally no speed increase, highlighting the 
sensitivity of this highly nonlinear effect. The Acb = 7r 
collisions remain elastic. For uj r ~ 2tt x 38 Hz, Nc de- 
creases to As and we cannot generate BSW solutions 
beyond this. If Ac is not exceeded all collisions should 
become elastic in the quasi-lD limit of large u r . 

In addition to radial trapping (u) r — 2?t x 17.5 Hz) we 
now add axial confinement lo z = 2?t x 6.84 Hz to compare 
with experiments Q. Note we have slightly modified u> z 
from the quoted experimental value of lj z — 2n x 6.8 Hz in 
order to match the experimental data, as we will see. It 
was observed that the number of atoms left in the conden- 
sate remnant was greater than Ac, and divided between 
two states, each with population As < A<c, which oscil- 
late axially. We continue to employ As = 1500, which 
satisfies these conditions. Note that the collapse-induced 
conversion of interaction into kinetic energy discussed 
above may explain the formation of oscillating BSWs in 
the experiment rather than stationary ones. The BSWs 
are positioned at zq = ±14 /im. We no longer apply 
a momentum kick since the trap accelerates the BSWs. 
The dynamics at early times for the cases of Acb = and 
Acb = it are shown in Fig. 0Ja) and (b) respectively Ini- 
tially the BSWs accelerate identically towards the trap 
centre. As they interact we see the different collisional 
dynamics, with the BSWs overlapping for Acb = and 
'bouncing' for Acb = %. However, at these early times the 
overall dynamics for both cases are similar: the BSWs 
emerge from the collisions with large collective excita- 
tions being generated. In Fig. 0Jc) we plot the early 
evolution of the axial full width half maximum FWHM 
(related to a gaussian fit) of the Acb = (dashed line) 
and Ac/> = tt (solid line) systems and compare with the 
experimental data Q (circles) . At early times there is lit- 
tle difference between the Acb = and 7r dynamics. The 
dominant feature in both the experimental and numerical 
data is an oscillation at 2u> z , corresponding to the centre- 
of-mass oscillations of the BSWs. We have chosen the ini- 
tial BSW separation to approximately match the maxi- 
mum experimental widths. The experimental widths de- 
crease to around 18/im which represents the resolution 
limit of the measurements. The numerical minima occur 
at FWHM ^2 — 5 (this is typically larger for Acp = tt 
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FIG. 4: (a)-(b) Evolution of the radially-integrated axial den- 
sity at early times for two BSWs (Ns = 1500 and initial am- 
plitude zo — ±14^tm) with (a) A(f> = and (b) A(f> = ir 
{lo z = 2tt x 6.84 Hz and ui r = 2tt x 17.5 Hz), (c) Full- width- 
half-maximum FWHM of a gaussian fit to the axial density for 
A(j> — (dashed line), A(j> — -n (solid line), and experimental 
data (points) 3]. Note we have time-shifted the experimen- 
tal data to match the phase of the oscillations at early times. 
(d)-(f) Same as (a)-(c) but for late times. 



due to the repulsive nature of the interaction). Whereas 
our simulations begin at maximum widths, the experi- 
mental data begins at a minimum 3] , indicating that the 
BSWs are created in close proximity. Note that in order 
to match the phase of the oscillations in Fig.^Jc) we have 
shifted the experimental data in time. 

We do not detect the secondary 2ui r frequency com- 
ponent observed experimentally. However, atoms ejected 
during the initial experimental collapse lead to a signifi- 
cant oscillation at 2io r Q ■ Furthermore we have verified 
that any offset of the BSWs from the z-axis will induce 
a 2cu r oscillation. The second largest frequency compo- 
nent in the simulations occurs at u> « 2tt x 27 Hz with an 
amplitude of only around 5% of the 2oj z component, and 
arises from the collective modes excited in the BSWs. 

Importantly, the experiment oscillations were observed 
for over 3 s with no discernable damping [3j . The dynam- 
ics at late times are shown in Fig.0Jd)-(f). The repeated 
collisions gradually destabilise the collisions through the 



growth of collective excitations. Partial collapse instabil- 
ities eject highly energetic atoms from the BSWs, which 
leads to the formation of a diffuse spread of low den- 
sity 'hot' atoms. The A(j> — case is most prone to 
instability, and by late times the BSWs have collapsed 
to a central highly-unstable wavepacket (Fig.^fc)). The 
FWHM for this case lies well below the experimental val- 
ues (Fig. ^Jf), dashed line). For A<fi = 7r, two oscillating 
wavepackets persist at late times, albeit in highly-excited 
states (Fig.ite)). The FWHM for A<j) = tt (Fig. Hp), 
solid line) remains close to the experimental values even 
at these late times. Note that in order to match the 
phase of the experimental data at late times we find it 
necessary to modify the quoted experimental axial trap 
frequency by around 0.5% to uj z — 2ir x 6.84 Hz. This 
agrees well with the fitted frequency in the experiment 
J3j, and implies that the BSW oscillations are an accu- 
rate measurement of trap oscillation frequency. 

In conclusion we find that bright solitary waves exhibit 
rich and complex behaviour, not present for ID solitons. 
We observe not just the upper critical population famil- 
iar with trapped attractive BECs but also a lower criti- 
cal population, below which the wavepacket is no longer 
capable of axially self-trapping. High-density collisions 
between BSWs induce instabilities. For large instability 
the BSWs are destroyed by a catastrophic collapse. For 
intermediate instability the collisions are inelastic, with 
higher or lower outgoing speed. We highlight the sensi- 
tive dependence of the instability on the phase difference 
between the BSWs, collisional speed, atomic population, 
and trap geometry, all of which are controllable exper- 
imentally. In particular elastic collisions are promoted 
for a 7r-phase difference between the BSWs and high col- 
lisional speeds. Furthermore, we find that the experi- 
mental observations of long-lived 'soliton' oscillations by 
Cornish et al. Q are successfully modelled in terms of 
two oscillating BSWs with a 7r-phase difference. 
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